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' We have investigated how the relative elemental abundances inferred from 

the solar upper atmosphere are affected by uncertainties in the dielectronic re- 
combination (DR) rate coefficients used to analyze the spectra. We find that the 
inferred relative abundances can be up to a factor of ~ 5 smaller or 1.6 times 
larger than those inferred using the currently recommended DR rate coefficients. 
We have also found a plausible set of variations to the DR rate coefficients which 
improve the inferred (and expected) isothermal nature of solar coronal observa- 
tions at heights of > 50 arcsec off the solar limb. Our results can be used to help 
prioritize the enormous amount of DR data needed for modeling solar and stellar 
upper atmospheres. Based on the work here, our list of needed rate coefficients 
for DR onto specific isoelectronic sequences reads, in decreasing order of impor- 
tance, as follows: 0-like, C-like, Be-like, N-like, B-like, F-like, Li-like, He-like, 
and Ne-like. It is our hope that this work will help to motivate and prioritize 
future experimental and theoretical studies of DR. 



Subject headings: atomic data - atomic processes - stars: abundances - stars: 
coronae - Sun: abundances - Sun: corona - Sun: transition region 
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1. Introduction 

It is now generally accepted that the elemental composition of the solar wind is different 
from that of the photosphere. Three decades of research have shown that in the slow speed 
solar wind (which is believed to originate from quiet coronal regions) the abundance ratio of 
low first ionization potential (FIP) elements (FIP < 10 eV) relative to higher FIP elements 
is larger than it is in the photosphere. The observed FIP enhancement is roughly a factor of 
four. With the advent of the Ulysses spacecraft in its polar orbit, it has become clear that 
in the fast solar wind (which emanates from coronal holes) the enhancement of the low FIP 
elements is certainly less than two, and is perhaps consistent with the solar photospheric 
abundance pattern (von Steiger, Geiss, & Gloeckler 1997; Feldman & Laming 2000). 

Similar phenomena are observed in spectroscopic studies of the solar corona and tran- 
sition region, which we will refer to as the solar upper atmosphere. Mapping the FIP en- 
hancement factor from the solar wind, through the corona and transition region, and into the 
photosphere (where by definition the factor is 1) is an important area of research for under- 
standing the formation of the solar wind. FIP factors for the solar wind are determined from 
in situ particle measurements and believed to be robust. FIP factors for the solar corona and 
transition region are inferred from spectroscopic observations and are sensitive to a number 
of uncertainties. Here we will explore how uncertainties in the dielectronic recombination 
(DR) rate coefficients limit our ability to infer FIP factors in the solar upper atmosphere. 

Of the many ionization and recombination rate coefficients which go into ionization 
balance calculations for solar and stellar coronae (i.e., electron- ionized plasmas), the high 
temperature DR rate coefficients are believed to be the most uncertain (Arnaud & Raymond 
1992; Mazzotta et al. 1998). Unlike most other atomic processes (e.g., ionization and excita- 
tion) where the direct contribution dominates the process, for DR is solely a resonant process. 
These resonances are doubly-excited, intermediate states which are highly correlated which 
makes the calculations theoretically and computationally challenging. 

Especially important are the high temperature DR rate coefficients for ions with par- 
tially filled L-shells. The astrophysical implications of these uncertainties are poorly under- 
stood. Investigating these implications is important because DR is the dominant electron-ion 
recombination process for most ions in cosmic plasmas. 

Typically, inferences of the FIP effect from the solar corona are made from spectra that 
exhibit lines from several charge states of the same element (see e.g., Malinovsky & Heroux 
1973; Laming et al. 1995; White et al. 2000). In such cases any uncertainties in the DR rate 
coefficients and the ionization balance calculations essentially cancel out. This is because 
an error that might increase the fraction of ions in a particular charge state, does so at the 
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expense of ionic fractions in neighboring charge states. If one observes a whole series of 
charge states to evaluate the FIP effect (e.g., Fe IX, X, XI, XII, XIII, XIV, XV, and XVI in 
the case of White et al. 2000), these problems are to a large extent obviated. 

Similar procedures are difficult to follow at solar transition region temperatures where 
fewer elements have substantial series of observable charge states. The situation is com- 
pounded by the temperature dependence of the emission measure, as discussed in § 3.2, and 
by the trend in solar physics instrumentation towards spectrometers with ever more limited 
bandpass while spatial and spectral resolution and sensitivity improve. For example, the 
spectrometers on the Solar Orbiting Heliospheric Observatory (SOHO) such as the Coronal 
Diagnostic Spectrometer (CDS; Harrison et al. 1995, 1997) and the Solar Ultraviolet Mea- 
surements of Emitted Radiation (SUMER; Wilhelm et al. 1995, 1997; Lemaire et al. 1997) 
must scan their detectors across the bandpass to build up a full spectrum. This limits their 
utility in observing time varying plasmas. The Extreme Ultraviolet Imaging Spectrometer 
(Culhane et al. 2001) will use multilayer coated gratings, which will further restrict the 
bandpass. The proposed wavelength ranges, ~ 170 — 210 and ~ 250 — 290 A will provide 
good coverage in strong lines of Fe IX - XVI and Fe XXIV, but are less satisfactory for ions of 
other elements. For these reasons accurate DR rate coefficients for the few lines that will be 
studied in current and future solar satellite missions have become more essential than ever. 

Here we investigate the extent to which FIP enhancement factors, inferred using indi- 
vidual line ratios, can be affected by uncertainties in the DR rate coefficients used in the 
data analysis. We hope to motivate further experimental and theoretical studies of DR. We 
begin by reviewing in § 2 the status of the theoretical DR results currently used in ionization 
balance calculations. We focus primarily on the rate coefficients for DR onto the L-shell ions 
Ne VI through Ne VIII, Mg VI through Mg IX, Si VI through Si XII, and S IX and S X. Lines 
from these ions are commonly observed from solar and stellar coronae. In § 3 we discuss how 
relative abundances are determined for solar and stellar upper atmospheres and the current 
observational status of relative abundance determinations. We describe in § 4 the ionization 
balance calculations, how the calculations are sensitive to our estimated uncertainties in the 
DR rate coefficients, and present some implications. We conclude in § 5 by making some 
specific recommendation for future directions in theoretical and experimental studies of DR. 



2. Dielectronic Recombination (DR) 

DR is a two-step process. It begins when a free electron collisionally excites an ion 
and is simultaneously captured (i.e., dielectronic capture). Core excitations of the ion can 
be labeled Nlj — > N'l'f. Here N is the principal quantum number of the core electron, 
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/ is its angular momentum, and j the total angular momentum of the core electron. In 
general, the most important DR channels are those via AN = N' — N = or 1 core 
excitations. The incident electron is captured into some nl" Rydberg level forming any one 
of an infinite number of intermediate doubly-excited states of the recombining ion. These 
states can either autoionize or radiatively stabilize, the latter of which completes the DR 
process. Energy conservation requires E\. = AE — E b , where Ek is the kinetic energy of 
the incident electron, AE the excitation energy of the core electron, and E b the binding 
energy released when the free electron is captured. AE and E b are quantized, making DR is 
a resonant process. The DR rate coefficient represents the convolution of these resonances 
with a Maxwellian electron energy distribution. 

Here we are interested in low-electron-density, zero-field DR rate coefficients. Ionization 
balance calculations generally assume that electron densities are low enough that collisional 
ionization of the weakly bound electron captured in the DR process is unimportant. These 
calculations also assume that any external electric and magnetic fields are too weak to affect 
the DR rate coefficient, and that metastable populations in all ions are insignificant. 

2.1. Overview of Theory and Experiment 

Reliable calculations of DR are extremely challenging theoretically and computationally. 
In contrast to most other atomic processes (e.g., ionization and excitation) where the direct 
contribution dominates the process, DR is solely a resonant process. These resonances 
are doubly-excited, intermediate states which are highly correlated. These states greatly 
complicate calculations, since an accurate treatment of electron correlation is required (Zong 
et al. 1997). Calculations also require accounting for an infinite number of states. This is 
clearly impossible using a finite basis expansion. Approximations must be made to make the 
calculations tractable (Hahn 1993). 

Many different theoretical techniques have been used to calculate DR rate coefficients 
for plasma modeling. In the past, semi-empirical expressions such as the Burgess (1965) for- 
mula along with modified versions by Burgess & Tworkowski (1976) and Merts et al. (1976) 
were derived to calculate DR rate coefficients. More recently, a number of sophisticated the- 
oretical approaches have been developed, among them configuration-averaging (Griffin et al. 
1985); single-configuration LS'-coupling (McLaughlin & Hahn 1984); quantum defect (Bell 
& Seaton 1985); intermediate-coupling (Badnell & Pindzola 1989b); nonrelativistic, multi- 
configuration Hartree-Fock (MCHF; Nilsen 1986; Schippers et al. 1998); semi-relativistic, 
multiconfiguration Breit-Pauli (MCBP; Badnell 1986); fully relativistic, multiconfiguration 
Dirac-Fock (MCDF; Chen 1985); and relativistic many-body perturbation theory (RMBPT; 
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Zong et al. 1997) methods. Other methods include the Hebrew University Lawrence Liv- 
ermore Atomic Code (HULLAC) which uses a relativistic, multiconfiguration, parametric 
potential method (Mitnik et al. 1998) and i?-matrix methods where radiative recombination 
(RR) and DR are treated in a unified manner in the close-coupling approximation (Nahar & 
Pradhan 1994; Robicheaux et al. 1995). However, the result of all these different theoretical 
techniques and required approximations are DR rate coefficients which often differ by factors 
of ~ 2 to 4 or more (Arnaud & Raymond 1992; Savin et al. 1997, 1999; Savin 1999, 2000). 

Cosmic plasmas are most commonly modeled using the recommended DR rate coef- 
ficients of Aldrovandi & Pequignot (1973), Shull & van Steenberg (1982), Nussbaumer & 
Storey (1983, 1984, 1986, 1987), Arnaud & Rothenflug (1985), Landini & Monsignori Fossi 
(1991), Arnaud & Raymond (1992), and Mazzotta et al. (1998). Essentially none of these 
DR rate coefficients have been calculated using state-of-the-art techniques (i.e. MCHF, 
MCBP, MCDF, RMBPT, HULLAC, or i?-matrix techniques which include the spin-orbit 
interaction). The vast majority of these have been calculated using single-configuration pure 
LS'-coupling, semi-empirical formulae, or isoelectronic interpolations. 

It is important to carry out calculations using techniques more sophisticated than LS- 
coupling because LS'-coupling calculations are known not to include all possible autoioniza- 
tion levels contributing to the DR process, due to not including the spin-orbit interaction. 
As a result, such calculations provide only a lower limit for the DR rate coefficient (Badnell 
1988; Gorczyca & Badnell 1996). For example, for lithiumlike ions, Griffin et al. (1985) 
and Belie & Pradhan (1987) discussed how LS'-coupling accounts for only two-thirds of all 
possible AN = recombining channels. Intermediate coupling calculations yield a DR rate 
coefficient for lithiumlike C IV which is 50% larger than the LS-coupling result. Recent ion 
storage ring measurements and relativistic many-body perturbation calculations have veri- 
fied the importance of these LS'-forbidden autoionizing resonances (Mannervik et al. 1998). 
As for the various semi-empirical formulae, Savin (1999) showed that a priori it is not pos- 
sible to know which of the semi-empirical formulae will yield a result close to the true DR 
rate coefficient and which will be off by a factor of 2. 

Due to the complexity of the theoretical descriptions for the DR process, it is almost 
impossible to say a priori which approximations in the calculations are justified and which 
are not. Laboratory measurements are needed to ensure that even state-of-the-art techniques 
produce reliable results. For example, initial Fe XVI MCBP results were a factor of ~ 2 larger 
than ion storage ring measurements (Linkemann et al. 1995). This discrepancy was later 
resolved (Gorczyca & Badnell 1996); but without laboratory measurements, the error in the 
theory would probably have gone undiscovered for many years. 

A related difficulty is that, for many of the ions and DR resonances important in electron- 
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ionized cosmic plasmas, there are very few experimental techniques capable of carrying 
out measurements at the accuracy needed to provide benchmarks useful for the theorists. 
Tokamak and theta pinch measurements suffer from factor of 2 uncertainties (Griem 1988). 
Crossed electron-ion beams techniques (Miiller et al. 1987; Savin et al. 1996) have been 
limited in the energy range accessible as well as by low signal rates. Single-pass, merged 
electron-ion beam techniques (Dittner et al. 1987; Andersen et al. 1990, 1992; Schennach et 
al. 1994) can cover a wider energy range but also have low signal rates. Due to a combination 
of poor statistics and poorly controlled external electric and magnetic fields, it is difficult 
to use most of these crossed- and merged-beams measurements to benchmark zero-field DR 
calculations. Using either technique it is also difficult to determine the ion beam metastable 
population. Hence these techniques cannot be used to carry out absolute measurements for 
many ions with partially filled shells (Badnell et al. 1991). 

The two state-of-the-art experimental techniques for studying DR are electron beam 
ion traps (EBITs) and heavy-ion storage rings with a merged electron-ion beam interaction 
region. Both techniques store the ions long enough for all metastable states to decay radia- 
tively to the ground state. Storage ring techniques measure absolute DR resonance strengths. 
EBITs produce relative resonance strength measurements which can then be normalized to 
RR or electron impact excitation theory. 

Until recently EBIT measurements were essentially limited to studying only closed shell 
systems (Beiersdorfer et al. 1992; Smith et al. 2000). This has changed with the recent 
studies of DR onto Fe XXI - Fe XXIV (Gu et al. 1999; Gu 2000). EBITs, however, are 
still not capable of reliably studying DR for collision energies < 670 eV (Wargelin et al. 
2001). These are the dominant DR channels for most ions of second and third row elements, 
especially for ions with partially filled shells. Work is underway to overcome this limitation. 

The majority of the storage ring measurements for astrophysically important ions have 
been carried out using TSR and CRYRING (Miiller 1995). Measurements of DR resonance 
strengths and energies have been reported using TSR from to 2200 eV (Kenntner et al. 
1995; Linkemann et al. 1995) and using CRYRING from to « 50 eV (DeWitt et al. 1995, 
1996). Storage rings are the optimal technique for studying DR resonances at energies not 
accessible to EBITs. This makes storage rings particularly well suited for studying DR for 
ions of second and third row elements. The external fields in the merged beams section are 
also extremely well controlled which allows reliable, essentially zero-field DR measurements 
to be carried out (Savin et al. 2000). One limitation, though, is that ions with a charge- 
to-mass ratios q/m < 0.1 cannot be accelerated in CRYRING and TSR to a velocity where 
electron capture from the rest gas in the ring becomes negligible. The resulting poor signal- 
to-noise ratio prevents measurements of DR for such ions (Schippers et al. 1998; Wolf 1999). 
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To date, EBIT and storage ring DR measurements have been limited in the isoelectronic 
sequences studied and the energy ranges covered. Few data sets exist which are comprehen- 
sive enough to benchmark high temperature DR theory. And for those ions studied over 
a wide enough energy range, the benchmarked theory has not then been used to calculate 
the DR rate coefficients for all astrophysically important ions isoelectronic to the measured 
ion. As a result, only a very small fraction of the DR rate coefficients used for modeling 
electron- ionized cosmic plasmas have come from benchmarked state-of-the-art theory. 

2.2. Uncertainties in the Recommended DR Rate Coefficients 

Here we are interested in how uncertainties in the published DR rate coefficients affect 
the ionization balance calculations used to analyze solar and stellar spectra. The ions of 
interest are listed in § 1. However, the calculated abundance for a given ion is affected by 
the atomic rate coefficients for the ions one lower and higher in charge state (and less so by 
ions two or more lower and higher). To account for this, we extend the list of ions of interest 
by one charge state on both the low and high end for each element considered. Below we 
review and estimate uncertainties in the published high temperature rate coefficients for DR 
onto Ne V - Ne IX, Mg V - Mg X, Si V - Si XIII, and S VIII - S XL 

In the absence of measurements, the best way to estimate the uncertainty in the theo- 
retical DR rate coefficients is to compare calculations along an isoelectronic sequence. This 
yields a conservative estimate. Savin et al. (1999) showed that published rate coefficients 
do not necessarily give reliable upper and lower limits for the range in which the true DR 
rate coefficient lies. The true DR rate coefficient may lie outside of these limits. Hence our 
approach here may, in fact, underestimate the uncertainties in the DR rate coefficients 

We compare the range of calculated DR rate coefficients to the recommended DR rate 
coefficients from Mazzotta et al. (1998) in order to estimate the uncertainties in the "state- 
of-the-art" for ionization balance calculations. Their DR rate coefficients are meant for use 
at temperatures between 10 4 and 10 9 K. Wherever possible we also use laboratory measure- 
ments to benchmark theory. Because experimental data often do not exist for the ions we 
are interested, we take into account measurements carried out on those isoelectronic ions 
which are closest in atomic number to the ions of interest. We also only take into account 
measurements of those resonances which are relevant to the DR rate coefficients of inter- 
est here. In our discussion below, the predicted temperature of formation for an ion in an 
electron-ionized plasma is taken from the results of Mazzotta et al. Also, estimated uncer- 
tainties for DR rate coefficients onto a given ion are quoted for temperatures near where the 
ion peaks in fractional abundance in an electron-ionized plasma. 
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Based on the comparisons described below, for the ions of interest we have put together 
a list of factors by which to scale the recommended DR rate coefficients of Mazzotta et al. 
(1998). These scale factors are listed in Table 1. The 2 6 x 3 3 possible combinations yield a 
total of 1728 sets of DR rate coefficients which we will use in § 4. 

2.2.1. Onto Heliumlike Ne IX and Si XIII 

LS-coupling calculations were carried out by Jacobs et al. (1977b), for only Si XIII, 
and by Romanik (1988). MCHF calculations were published by Nilsen (1986) and Karim 
& Bhalla (1989). Chen (1986) has published MCDF rate coefficients. The results of Chen, 
Nilsen, Romanik, and Karim & Bhalla are in excellent agreement while the results of Jacob 
et al. are a factor of ~ 2 smaller. Mazzotta et al. (1998) use the rate coefficients of Chen. 

EBIT and electron beam ion source studies of DR have been carried out for a number 
of heliumlike ions. The most relevant measurements for our comparison here have been on 
Ne IX (Wargelin, Kahn, & Beiersdorfer 2000) and Ar XVII (Ali et al. 1991; Smith et al. 1996). 
There is an estimated total experimental uncertainty of m 20% for these results. Overall 
there is good agreement between theory and experiment. Calculations using the technique 
of Karim & Bhalla (1989) are in good agreement with the Ar XVII measurements of Ali et 
al. (1991). Calculations using the technique of Chen (1986) are in good agreement with the 
results on Ar XVII (Smith et al. 1996). For Ne IX, for the sum of the measured Ne IX KLL 
resonances, experiment lies a factor of 1.16 below calculations using the technique of Karim 
& Bhalla (1989), a factor of 1.23 below calculations using the technique of Chen (1986), and 
a factor of 1.32 below calculations using the technique of Nilsen (1986). 

Overall the calculations of Chen (1986) agree with laboratory results to within m 20%. 
We take this estimate as the uncertainty in the recommended rate coefficients of Mazzotta 
et al. (1998) for DR onto Ne IX and Si XIII. 

2.2.2. Onto Lithiumlike Ne VIII, Mg X, and Si XII 

LS'-coupling DR calculations have been published by several different workers. Jacobs 
et al. (1977b, 1979) and Romanik (1988) published coefficients for Ne VIII, Mg X, and Si XII. 
Roszman (1987a) gave rate coefficients for Ne VIII. MCDF results were published by Chen 
(1991) for Ne VIII and Si XII. For reasons discussed in § 2.1, we do not use the published 
LS results to estimate the uncertainty in the DR rate coefficients onto lithiumlike ions. 

Measurements have been carried out for a number of lithiumlike ions. Storage ring 
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measurements have been carried out on Ne VIII (Zong et al. 1998) and Si XII (Kenntner 
1995; Kenntner et al. 1995; Bartsch et al. 1997). Uncertainties were typically pa 20%. 

At the temperature of peak formation, the dominant DR channel for lithiumlike ions up 
to Si XII is via AN = core excitations. For Si XII the AN = and AN = 1 (N = 2 -> N' = 
3) contributions are comparable (Chen 1986). For AN = DR onto Ne VIII MCBP theory 
was ~ 20% below the experimental data. MCPB theory was ~ 10 — 20% below experiment 
for AN = DR onto Si XII. For Si XII AN = 1 DR via iV = 2 -> N' = 3 core excitations, 
MCBP theory was larger than experiment by pa 10 — 20% for the ls 2 3/3/' resonances. For 
the ls 2 3/n/' (n > 4) resonances in Si XII, theory lies slightly below experiment. These 
comparisons suggest that the accuracy of the MCBP results is ~ ±20%. 

None of the measured resonance strengths have been compared with results from the 
theoretical techniques used to calculate the published rate coefficients for Ne VIII, Mg X, 
and Si XII. This makes it difficult to use the measurements to infer an uncertainty in the 
recommended DR rate coefficients. We try to do this indirectly by comparing MCBP results 
with published MCDF results. The MCDF results of Chen (1991) agree to within « 10% 
with the MCBP rate coefficients of Badnell & Pindzola (1989b) for O VI and of Badnell & 
O'Mullane (1999) for Ar XVI. Considering the comparisons between the MCBP and MCDF 
results and given the estimated uncertainty in the MCBP results, we therefore estimate that 
there may be an pa ±20% uncertainty in the MCDF rate coefficients of Chen (1991). 

Mazzotta et al. (1998) use the rate coefficients of Chen (1991) and interpolate isoelec- 
tronically for those ions which Chen did not calculate. Chen did not present results for 
T < 10 5 K. As a result the fitted rate coefficients of Mazzotta et al. do not have the correct 
low temperature behavior. This can readily be seen by plotting their recommended rate 
coefficient for C IV and comparing it with published C IV rate coefficients (Schippers et 
al. 2001). This error in the rate coefficients of Mazzotta et al. will affect ionization bal- 
ance calculations for photoionized plasmas but is expected to have little effect on modeling 
electron-ionized plasmas. Given the estimated uncertainty in the results of Chen, we esti- 
mate that for Ne VIII, Mg X, and Si XII there is an pa ±20% uncertainty in the relevant DR 
rate coefficients of Mazzotta et al. 



2.2.3. Onto Berylliumlike Ne VII, Mg IX, and Si XI 

LS'-coupling rate coefficients were published for Ne VII Mg IX, and Si XI by Jacobs et 
al. (1977b, 1979) and Romanik (1988). Near the temperatures of peak formation, their high 
temperature results are in good agreement. MCBP results were presented by Badnell (1987) 



-10- 



for Ne VII, Mg IX, and Si XI. At the temperature of peak formation, the rate coefficients of 
Jacobs et al. and Romanik are a factor of rs 1.4 — 1.6 times larger than those of Badnell. 

Mazzotta et al. (1998) use the rate coefficients of Badnell (1987) and appear to interpo- 
late along the isoelectronic sequence for those ions which Badnell did not calculate. Badnell 
does not present any results for T < 10 5 K. A comparison of the rate coefficients from Maz- 
zotta et al. with those from Romanik (1988) shows that the former do not have the correct 
low temperature behavior. This will be important in modeling the ionization structure of 
photoionized plasmas but is not an issue here. 

In the absence of laboratory benchmarks, we use the various published DR calculations 
to provide upper and lower limits for the DR rate coefficients. As discussed above, this is a 
conservative estimate of the uncertainty. We estimate that for Ne VII, Mg IX, and Si XI, the 
uncertainty in the relevant DR rate coefficients of Mazzotta et al. (1998) is +60% and -0%. 

2.2.4. Onto Boronlike Ne VI, Mg VIII, and Si X 

Jacobs et al. (1977b, 1979) reported LS-coupling DR rate coefficients for Ne VI, Mg VIII, 
and Si X. Nahar (1995) used R- matrix techniques to calculate RR+DR rate coefficients in 
LS'-coupling for Ne VI, Mg VIII, and Si X. Mazzotta et al. (1998) recommends DR rate 
coefficients based on the calculations of Nahar after apparently subtracting the theoretical 
RR rate coefficient from her results. 

Given the paucity of theoretical calculations and appropriate laboratory measurements, 
it is difficult to estimate the uncertainty in the theoretical DR rate coefficients for the ions 
of interest. We attempt to do this indirectly by comparing the results for DR onto the 
isoelectronic ions C II, N III, and O IV for which a number of different calculations exist. 
LS'-coupling rate coefficients were reported by Jacobs et al. (1978) and Ramadan & Hahn 
(1989). Fits to the results of Jacobs et al. (1978) were reported by Shull & van Steenberg 
(1982). i?-matrix results using LS-coupling were presented by Nahar (1995). Badnell & 
Pindzola (1989a) carried out MCBP calculations and Safronova & Kato (1998) carried out 
MCHF calculations. The various theoretical techniques used have not yet converged to 
the same rate coefficients. Mazzotta et al. (1998) recommends the rate coefficients of Nahar 
(1995). At the temperatures of peak formation for C II, N III, and O IV, the other theoretical 
rate coefficients are a factor of ~ 1.0 — 1.7 larger than the recommended rate coefficients. 
We estimate that for Ne VI, Mg VIII, and Si X, the uncertainty in the relevant DR rate 
coefficients of Mazzotta et al. is +70% and -0%. 
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2.2.5. Onto Carbonlike Ne V, Mg VII, Si IX, and S XI 

There have been very few calculations of high temperature DR for the ions of interest 
here. LS'-coupling results have been presented for Ne V, Mg VII, Si IX, and S XI by Jacobs 
et al. (1977b, 1979). Mazzotta et al. (1998) use the rate coefficients of Jacobs et al. scaled up 
by a factor of m 1.6. This factor was apparently derived by scaling the isoelectronic Fe XXI 
results of Jacobs et al. (1977a), calculated using the same technique as Jacobs et al. (1977b, 
1979), to the recommended DR rate coefficient of Arnaud & Raymond (1992). 

Considering how few relevant calculations exist, it is difficult to estimate the uncertainty 
in the DR rate coefficients for the ions of interest here. We attempt to do this indirectly 
using theoretical results for high temperature DR onto N II and O III. Jacobs et al. (1978) 
presented LS-coupling results for N II and O III which were fitted by Shull & van Steenberg 
(1982). i?-matrix results in LS'-coupling were calculated by Nahar & Pradhan (1997) for N II 
and by Nahar (1999) for O III. MCBP and LS rate coefficients were calculated by Badnell & 
Pindzola (1989b) and Roszman (1989), respectively, for O III. Note that there is an apparent 
error in Table 1 of Roszman (1989). We had to reduce the DR rate coefficients by an order of 
magnitude in order for the tabulated rate coefficients to match those shown in Figures 1 and 
2 of his paper. Near the high temperature DR peak for N II, the rate coefficient of Nahar 
& Pradhan (1997) is a factor of « 2.2 smaller than that of Jaocbs et al. For O III near this 
peak, the rate coefficient of Roszman is a factor of « 1.2 times larger than that of Jacobs 
et al. The rate coefficients of Badnell & Pindzola and Nahar & Pradhan are, respectively, 
~ 1.25 and « 1.7 smaller than the results of Jacobs et al. 

For DR onto Ne V, Mg VII, Si IX, and S XI, we use the recommended rate coefficients 
of Mazzotta et al. (1998). A reduction in these rate coefficients by 38% brings their results 
into agreement with the original results of Jacobs et al. Considering the comparison of the 
various theoretical DR rate coefficients of N II and O III, a reduction in the rate coefficients 
of Mazzotta et al. by 69% represents our estimated lower limit to the uncertainty in their 
recommended DR rate coefficients. 



2.2.6. Onto Nitrogenlike Mg VI, Si VIII, and S X 

We are aware only of the LS'-coupling calculations by Jacobs et al. (1977b, 1979) for 
Mg VI, Si VIII, and S X. Mazzotta et al. (1998) use the rate coefficients of Jacobs et al. For 
S X, Mazzotta et al. do not include the low temperature results of Jacobs et al. This is not 
an issue here but may be important in photoionized plasmas. 

The paucity of calculations for DR onto nitrogenlike ions makes it difficult to estimate 
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the uncertainty in the published theoretical rate coefficients. We attempt to do this indirectly 
using published rate coefficients for DR onto the isoelectronic O II. High temperature DR 
rate coefficients have been calculated using LS'-coupling by Jacobs et al. (1978) and Terao 
et al. (1991). Fits to the results of Jacobs et al. were published by Shull & van Steenberg 
(1982). Nahar (1999) published i?-matrix rate coefficients (e.g. RR+DR) using LS'-coupling. 
Intermediate coupling calculations have been given by Badnell & Pindzola (1989b) which 
were later improved upon by Badnell (1992). 

The results of Jacobs et al. (1978) and Badnell & Pindzola (1989b) are in good agree- 
ment. However, Badnell (1992) recalculated the DR rate coefficient taking into account 
correlation between the n = 2 and n = 3 shells. His new rate coefficient is a factor of ~ 1.6 
times smaller than that of Jacobs et al. (1978) and Badnell & Pindzola (1989a). Terao et 
al. (1991) revised their published rate coefficients upward by « 20% (Badnell 1992). This 
revised result is a factor of 2.2 smaller than the results of Jacobs et al. (1978). To de- 
termine the DR rate coefficient of Nahar (1999), we subtracted out the estimated RR rate 
coefficient. By extrapolating her recombination rate coefficient at low temperatures, where 
DR is unimportant, we estimated the RR rate coefficient at the relevant high temperatures. 
The resulting DR rate coefficient lies a factor of ~ 2.2 below that of Jacobs et al. Taking 
into account these comparisons, we estimate the uncertainty in the relevant rate coefficients 
of Mazzotta et al. for DR onto Mg VI, Si VIII, and S X to be +0% and -55%. 



2.2.7. Onto Oxygenlike Mg V, Si VII, and SIX 

Jacobs et al. (1977b, 1979) have published LS rate coefficients for Mg V, Si VII, and 
S IX. These rate coefficients were fitted by Shull & van Steenberg (1982). Mazzotta et al. 
(1998) use these results. We estimate the uncertainties in the DR rate coefficients for these 
ions by comparing the rate coefficients for DR onto the isoelectronic Fe XIX for which a 
number of different calculations exist. Jacobs et al. (1977a) and Roszman (1987c) published 
LS-coupling results. Jacobs et al. (1977a) used the same technique as Jacobs et al. (1977b, 
1979). Their results were fitted by Shull & van Steenberg (1982) for use in ionization balance 
calculations. Hartree-Fock calculations with relativistic corrections, using the code of Cowan 
(1981), were reported by Dasgupta & Whitney (1994). 

In making our comparisons between the different calculations, we need to take into 
account the relative importance of the different DR channels for the ions of interest. For 
Fe XIX, near the temperature of peak formation the AN = 1 (N = 2 — > N' = 3) channel 
appears to dominate the DR process (Roszman 1987c; Dasgupta & Whitney 1994). It is 
unclear how this extrapolates to the less highly charged Mg V, Si VII, and S IX. 
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For Fe XIX at the temperature of peak formation, the rate coefficients of Roszman 
(1987c) and Dasgupta & Whitney (1994) are, respectively, factors of approximately 2.7 and 
3.6 times larger than the result of Jacobs et al. (1977a). Mazzotta et al. use the rates of 
Jacobs et al. (1977b, 1979) for Mg V, Si VII, and S IX. Based on our comparison for Fe XIX, 
we estimate the uncertainty in the high temperature rate coefficients for DR onto these ions 
to be +260% and -0%. 



2.2.8. Onto Fluorinelike Si VI and S VIII 

Jacobs et al. (1977b, 1979) carried out LS'-coupling calculations for Si VI and S VIII. 
We are unaware of any other calculations for these ions. Shull & van Steenberg (1982) fit 
the data of Jacobs et al. Mazzotta et al. (1998) use these fits. 

We can estimate the important DR channels in Si VI and S VIII and the uncertainties 
in the relevant DR rate coefficients using the theoretical results for the isoelectronic ions 
Fe XVIII. Jacobs et al. (1977a), using the same techniques as Jacobs et al. (1977b, 1979), 
have published an LS rate coefficient. Shull & van Steenberg (1982) presented a fit to 
this rate coefficient for use in plasma modeling. Roszman (1987b) calculated an LS rate 
coefficient. Dasgupta & Whitney (1990) published Hartree-Fock calculations with relativistic 
corrections using the code of Cowan (1981). MCDF results were reported by Chen (1988a). 
At the temperatures of peak formation for Fe XVIII, the AiV = 1 (N = 2 — > N' = 3) channel 
dominates the DR process. It is likely that the same situation exists for Si VI and S VIII. At 
the temperature of peak formation for Fe XVIII, the rate coefficients of Roszman (1987b), 
Chen (1988a), and Dasgupta & Whitney (1990) lie factors of approximately 3.5, 4.2, and 
4.7, respectively, above the results of Jacobs et al. (1977a). Here, for Si VI and S VIII, we 
use the recommended high temperature DR rate coefficients of Mazzotta et al. (1998) and 
estimate the uncertainty to be +370% and -0%. 



2.2.9. Onto Neonlike Si V 

Jacobs et al. (1977b) and Romanik (1988) have calculated DR rate coefficients using 
LS'-coupling for Si V. At the peak in the DR rate coefficient, the results of Jacobs et al. lie a 
factor of ~ 1.9 below those of Romanik. Mazzotta et al. (1998) use the results of Romanik. 
We use the theoretical rate coefficients for DR onto Ar IX to estimate the uncertainty in 
the recommended Si V data. LS rate coefficients were given by Romanik. Chen (1986) 
published MCDF rate coefficients. HULLAC results were given by Fournier et al. (1997). 
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The results of Romanik (1988) and Fournier et al. (1997) are in excellent agreement. 
The calculations of Chen (1986) lie ~ 15% below those of Romanik and Fournier et al. The 
close agreement between these three different calculations and the large difference between 
the results of Jacobs et al. (1977b) and Romanik for Si V, strongly suggests an error in the 
reported results of Jacobs et al. This is supported by a comparison of rate coefficients for 
DR onto Fe XVII by Arnaud & Raymond (1992). They found that the results of Chen and 
Romanik were in good agreement but that the results of Jacobs et al. (1977a), calculated 
using the same techniques as Jacobs et al. (1977b), were a factor of ?a 4.7 times smaller. 
Considering these various comparisons, we use the recommended rate coefficient of Mazzotta 
et al. (1998) for Si V and estimate the uncertainty to be +0% and -15%. 

3. Solar and Stellar Upper Atmosphere Abundance Observations 



Spectral line ratios can be used to determine relative abundances in cosmic plasmas. 
Photons are emitted in a spectral line at a rate per second per steradian given by 



where R is the rate, A-^ the radiative decay rate from excited ionic level j to i, n- } the number 
density of level j, n q the number density for the g-times charged ion in question, the 
number density for the element from which this ion is formed, nn the hydrogen number 
density, n e the electron number density, and AV the emitting volume. The radiation is 
emitted into 4n steradians. For full disk spectra it is usual to include a factor of 1/2 to 
account for photons emitted towards the solar disk and absorbed by the photosphere. At 
low densities (i.e., in coronal plasma), we can write rij/n q = Cfi c /A^ where Cj is the total 
excitation rate coefficient of the excited level j by electron collisions and we have assumed 
the branching ratio for the j — > i radiative transition to be 1. For the transitions of interest 
in the solar FIP diagnostics under consideration here (i.e. n = 2 — > n = 2 transitions in 
L-shell ions), excitation rate coefficients due to other processes are negligible compared to 
those from electron impact excitation. Defining f q = n q /nA and changing the integration 
variable from volume to temperature gives 



where parameters not dependent on the electron temperature have been taken outside the 
integral. The quantity n 2 e dV/dT is the differential emission measure (DEM, also sometimes 



3.1. Determining Relative Abundances 




(1) 




(2) 
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defined as n 2 c dV/d log T). The usual procedure in determining the DEM is to assume some 
functional form (usually DEM oc T a ), evaluate R for all lines, compare these values with 
observations and iterate on the DEM distribution until satisfactory agreement is obtained. 
This can be done for all the lines of one particular element (usually Fe). Then lines from 
other elements can be analyzed the same way to determine relative element abundances. 
However in this procedure the role of uncertainties in particular atomic rate coefficients 
becomes obscured by the quantity of data and number of iterations required. 

A second technique for determining relative abundances is to use a line from an ion of 
a specific element and another line from an ion of a different element. This is the technique 
which will most likely be used to infer relative abundances from future solar satellite missions. 
The lines are chosen so that the respective values of Cj/ q , to a first approximation, have 
similar dependences on temperature, differing only by a multiplicative constant. Using 
Equation 2 to take the ratio of the two lines, it is easy to show that the ratio of the two 
integrals reduces to this constant, which is determined from atomic parameters, times the 
relative abundances. In this approximation, the shape of the DEM has no affect on the 
results. This technique has the advantage of allowing one to study how uncertainties in 
atomic physics affect the inferred relative abundances. 



3.2. Solar Disk Observations 

The first systematic analysis of coronal abundances over a wide temperature region using 
high quality spectral data was probably that of Laming et al. (1995) who used a variant of 
the first DEM method described above. They analyzed the full-disk solar corona spectrum 
of Malinovsky & Heroux (1973) and covered a range of logT = 5.5 — 6.5. They found an 
abundance pattern consistent with that of the slow speed solar wind, at least for solar coronal 
plasma with an electron temperature logT > 5.8 (i.e. T > 9 x 10 5 K). More recent studies 
with lines from ions formed at similar high temperatures corroborate this result (Laming et 
al. 1999; White et al. 2000). But in what was a surprise at the time, they found essentially 
no FIP enhancement in transition region plasma at logT < 5.8. A similar result was found 
earlier by Noci et al. (1988), who used data from the Skylab S-055 EUV spectroheliometer. 

Interpreting spectra from these cooler temperatures poses a number of problems. This 
may be in part due to complications in the solar physics at these temperatures (Feldman 
& Laming 1994), or due to problems in the interpretation of the data. At transition region 
temperatures, the solar emission measure EM = n 2 e V (T) is an increasing function of T. 
The "canonical" quiet sun behaviour is EM oc T 1,5 (Jordan 1980), though in certain solar 
features power laws of T 4 or higher may be appropriate (Cargill 1994; Cargill & Klimchuk 
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1997). This is quite different from the situation at logT > 6 where in the quiet sun the 
emission measure is in general much less steep, and may even be flat in certain temperature 
regimes (e.g., Laming et al. 1995). One consequence of the EM temperature dependence is 
that in observations of ions formed at logT > 5.8, where the emission measure is not sig- 
nificantly temperature dependent, the detected line emission comes primarily from regions 
near the temperatures of peak fractional abundance for the observed ions. However, this is 
not the case for ions with formation temperatures of logT < 5.8. For these ions, the steeply 
rising EM skews the temperature range over which emission from the ions is sampled to 
temperatures well above those where the ions peak in fractional abundance. Under these 
conditions uncertainties in the atomic data used to calculate the ionization balance can sig- 
nificantly affect the inferred properties of the observed plasma. In addition, at temperatures 
significantly above or below those where an ion is formed in ionization equilibrium, we find 
that the fractional ion abundances calculated by balancing the electron-ion ionization and 
recombination rate coefficients are extremely sensitive to errors in the atomic data. 

Much of the work prior to that of Laming et al. (1995) had clearly observed FIP effects 
at logT < 5.8, principally using the Mg VI/Ne VI and Ca IX/Ne VII intensity ratios in 
the second, simpler method outlined in § 3.1 for determining relative abundances (Feldman 
& Widing 1990, 1993; Widing & Feldman 1989, 1992, 1993). One reason for the apparent 
disagreement between these studies and that of Laming et al. (1995) is probably that while 
Laming et al. studied a full disk integrated spectrum, the prior work had concentrated on 
discrete solar features observed by the Skylab SO-82A spectroheliograph. Specifically, the 
Skylab investigations had concerned a coronal polar plume (Widing & Feldman 1992), an 
impulsive flare (Feldman & Widing 1990), an open field active region (Widing & Feldman 
1993), or a variety of such features (Widing & Feldman 1989). 

Interestingly the study of a coronal hole using a Mg VI/Ne VI line ratio (Feldman & 
Widing 1993) yielded an apparent FIP enhancement of about 2-2.5 at logT 5.6. This is 
slightly higher than the findings of Laming et al. (1995) for the full disk sun in this lower 
temperature range. However, given the estimated uncertainties of 50% in the Feldman & 
Widing (1993) and Laming et al. (1995) in this temperature range, and the different methods 
of analysis, there is a suggestion that these results are in agreement. The FIP enhancements 
of Feldman & Widing (1993) are also similar to or slightly higher than those observed in 
the fast speed solar wind, which is believed to emanate from coronal holes. The results of 
Feldman & Widing (1993) and Laming et al. (1995) support the hypothesis that the FIP 
factor decreases as one moves inward from the solar corona to the transition region. 

A considerable amount of attention has been given to the open field active region re- 
ported by Widing & Feldman (1993) and the coronal polar plume reported by Widing & 
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Feldman (1992). This is due to the observed FIP enhancements being greater than an order 
of magnitude. SOHO observations of open field active regions using CDS have found results 
(Young & Mason 1997) similar to those of the active region of Widing & Feldman (1993). 
SOHO observations have also been carried out of coronal polar plumes, principally using 
SUMER. Doschek et al. (1998) made observations of the Si/Ne abundance ratio (as a proxy 
for the FIP effect) and found a FIP enhancement of only a factor of 2. Other SUMER 
observers find similar results (del Zanna & Bromage 1999; Young et al. 1999). So far, it has 
not been possible to reproduce the polar plume results of Widing & Feldman (1992). 

The emission measure in the open field active region analyzed by Widing & Feldman 
(1993) is relatively flat with temperature and the inferred FIP enhancements of these features 
are believed to be reliable. However, the emission measure plotted by Widing & Feldman 
(1992) for their observed coronal polar plume rises very steeply with temperature. Although 
the Mg VI and Ne VI ionization fractions match each other very well near their temperature 
of peak formation, at higher temperatures the Mg VI fraction becomes larger than that for 
Ne VI. This is illustrated in Figure 1 where the relative intensities are plotted for the Mg VI 
1190.09 and Ne VI 558.59 A lines, using the ionization balance of Arnaud & Rothenflug 
(1985). The excitation rate coefficients are taken to vary as exp (—AE/kT) /VT, where 
AE is the excitation potential and T is the electron temperature. Under conditions where 
the temperature structure of the plasma weights the emission from these lines so that the 
temperature region greater than 5 x 10 5 is dominant, apparent abundance enhancements of 
Mg over Ne may result simply due to the relative ionization fractions. Realistic calculations 
by Doschek & Laming (2000) for the emission measure distributions found by Widing & 
Feldman (1992) suggest that an extra factor of 2-3 in the intensity ratio, on top of the usual 
FIP effect, could result from this effect. Thus if the true FIP effect is a factor of 2-4, an 
apparent abundance enhancement of an order of magnitude or more could result. 



3.3. Stellar Observations 

All observations of stellar coronae are necessarily disk integrated. In general in spectra 
of stellar coronae, individual line ratios for FIP effect diagnostics are not as readily available 
as they are in solar spectra. Most abundance results have been derived using these spectra 
from emission measure plots, and related techniques, where the effect of uncertainties in DR 
is much less transparent. However two papers that do use specific line ratios are Drake et 
al. (1995) on Procyon and Drake et al. (1997) on a Cen. These sources are perhaps the two 
best observed stellar coronae with the Extreme Ultraviolet Explorer (EUVE). 

In general, analyses on EUVE spectra of solar-like stars have found coronal abundance 



-18- 



anomalies similar to the solar FIP effect. Only Procyon appears to have a photospheric 
abundance corona. Coronal abundances in active binary stars do not show such effects, 
but hint at a trend of decreasing coronal metallicity with increasing activity. These and 
other results prior to the launch of the Chandra and XMM-Newton satellites are reviewed 
by Feldman & Laming (2000). With the advent of high resolution stellar x-ray spectroscopy 
with Chandra and XMM-Newton, the field of stellar coronae and their abundances is set to 
become a rich area of research. Already new results have raised the possibility of yet more 
varied coronal abundance patterns, e.g., the inverse FIP effect (Brinkmann et al. 2000) or 
enhanced Ne and Ar (Drake et al. 2001) in the coronae of HR 1099. 

4. Effects of DR Uncertainties on Relative Abundance Determinations 

Our work focuses on studying particular line ratios because of their importance for 
determining relative abundance ratios using data collected from past, current, and future 
solar satellite missions. These line ratios also offer a relatively transparent insight into 
the effect of uncertainties in the DR rate coefficients on abundance measurements. SOHO 
studies have been carried out using Mg VI/Ne VI, Mg VII/Ne VII, Mg IX/S IX, Mg IX/S X, 
Si IX/S IX, and Si IX/S X line ratios. Wavelengths for the observed lines are given in detail 
in Table 3. An additional abundance ratio that is of interest is Si X/S X, derived from the 
Si X 258.40 and 261.27 A lines and S X transitions at 259.50 and 264.24 A which has been 
used in EUVE abundance studies on the coronae of a Centauri (Drake et al. 1997). We 
consider all of the 7 above listed line ratios in our discussion below. 



where Ci ons , C rr>q , and Cd r , q are the rate coefficients for electron impact ionization, RR, and 
DR respectively, out of the charge state q. For the coronal plasmas under consideration 
here three-body recombination can be safely neglected. In ionization equilibrium, df q /dt = 
which gives Z independent linear equations for the Z + 1 charge states. The constraint that 
all charge state fractions must add up to unity supplies the final equation, allowing us to 
solve the Z + 1 linear equations by LU decomposition (see Press et al. 1992, for details) of 
the matrix formed by the right hand side of Equation 3. This set of linear equations is solved 
repeatedly with rate coefficients appropriate to different electron temperatures to find the 
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Ionization Balance Calculations 




(3) 
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temperature dependence of the different charge state fractions. Initially the ionization, RR 
adn DR rate coefficients are taken from the same sources as used by Mazzotta et al. (1998). 



4.2. Relative Abundances 

Using Equation 2, the relative abundances for two elements, Al and A2 is given by 

n Al = -RA1 x fAT C i 2 fq 2n elT dT 

n A2 RA2 LrC^n^dT 

This equation is first evaluated with ionization balance calculations which use the unsealed 
DR rate coefficients of Mazzotta et al. (1998). To investigate the effects that the various sets 
of DR rate coefficients have on the inferred value of n^i/n^ we define the quantity 

s= ( !&T C vf^ 2n l% dT \ / ( fAT C vft 2n e% dT \ ^ 

VJat Cjifq 1 ^ -jFpdT J ncw / \ f Ar Cji/A^^f dT / i d 

where the subscripts old and new refer, respectively, to ionization balance calculations carried 
out using the unsealed and scaled DR rate coefficients of Mazzotta et al. (1998). Multiplying 
the right-hand-side of Equation 4 by S yields the new inferred relative abundances. We note 
that a decrease (increase) in the calculated value of S corresponds to a decrease(increase) in 
the inferred FIP enhancement for a given value of Rai/Ra2- 



4.3. Effects 

For our calculations of S we made a number of simplifying assumptions. For the pairs 
of lines we are concerned with here, to a first approximation the values of Cj/ q have the 
same temperature dependence. To simplify the calculations we have assumed a flat DEM 
(oc T°). We further simplify the calculations by neglecting the temperature dependence of 
Cj, setting it equal to 1 for each ion. In this way we are able to focus specifically on how the 
uncertainty in the DR rate coefficients affects S. We then calculated S for all 1728 sets of 
DR rate coefficients. The calculations were carried out over the temperature range AT from 
logT = 5.0 — 7.0. This covers the temperature range over which the ions of interest form. For 
brevity, we have listed in Table 2 only the range of S values for all 1728 variations. Clearly 
the uncertainties in the DR rate coefficients can have a dramatic affect on any inferred FIP 
enhancement factors. We note that Table 2 does not display the correlations between the 
values of S for different line ratios for each set of DR variations. 
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We have reduced the set of DR variations using our current hypothesis for the structure 
of the FIP effect. Mg VII and Ne VII are transition region temperature ions which have been 
observed by Laming et al. (1999) from close to the solar limb out into the corona (where 
these ions are observed far away in temperature from where they peak in abundance). The 
inferred FIP effect is expected to be ~ 4 in the coronal observations and decrease as one 
moves to transition region temperatures closer to the solar limb. This expected behavior 
for the Mg/Ne abundance ration can be seen in Table 3 which is reproduced here from 
Laming et al. (1999). Line emission from Mg IX and S X is always dominated by coronal 
temperature plasmas and the Mg IX/S X ratio is expected to display a FIP factor of 4 at 
all positions listed in Table 3. Laming et al. (1999) measured a FIP factor of 3. For 
our selection criteria, we take only those DR variations for which the FIP effect inferred 
using Mg VII/Ne VII line ratio varies by less than ±25% and which also increase the FIP 
effect inferred using Mg IX/S X. The 25% is based on the estimated errors in the measured 
line ratios. We found 548 sets which met these constraints. This set was reduced to 274 
by taking into account the fact that changes in the DR rate coefficient onto neonlike Si V 
had essentially no affect on the results. We have listed in Table 2 the range of S for this 
reduced set of variations. Clearly, even in this reduced set the uncertainties in the DR rate 
coefficients can still have a dramatic affect on many inferred FIP enhancement factors. 



4.4. Solar Off-Limb Observations: A Further Test 



At heights of 50 arcsec or more from the solar limb, the corona becomes essentially 
isothermal. As one away from the solar limb, the amount of cooler transition region plasma 
decreases quickly with height, while the hotter coronal gas diminishes much less rapidly, 
and quickly becomes the dominant component. A striking illustration of this in the case 
of a solar equatorial streamer is given in Feldman et al. (1999). In Figure 3 of their paper 
they plot the variation of intensity with distance from the solar limb for various transitions 
from ions of Si VH-Si XII. All line intensities / show the same slope dlogl/dr and hence 
the same temperature, which from the slope evaluates to about 1.5 x 10 6 K. Feldman et al. 
plot in their Figure 4 the loci of emission measures determined from these Si ions. Possible 
calibration uncertainties exist for the Si XII results as it is at the extreme short wavelength 
end of the SUMER second order bandpass (Laming et al. 1997). Ignoring the Si XII loci, 
the temperature determined from the intersection of the remaining emission measure loci is 
logT = 6.11 ± 0.04. This is similar to that determined from the height dependence of the 
emission in the various lines, assumed to fall off according to hydrostatic equilibrium. 
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Feldman et al. (1999) used the older ionization balance of Arnaud & Rothenflug (1985). 
In order to assess the impact of our suggested changes to dielectronic recombination rate 
coefficients, we have remeasured the line intensities and replotted the emission measure 
determined using the ionization balance of Mazzotta et al. (1998). These results are shown in 
Figure 2. Comparing with Figure 4 of Feldman et al. (1999) we can see that the temperature 
at which the various curves intersect is now given by logT = 6.14 ± 0.05. 

We have investigated the effect of varying the DR rate coefficients in the Si ionization 
balance calculations. We have looked at variations of the rate coefficients in a "trial and 
error" fashion. Nine plausible sets of variations are given in Table 4. Others exist among our 
1728 attempts, but these nine serve to illustrate one effect we are interested in. In Figure 3(a) 
we replot the temperature region 6.0 < logT < 6.3 from the previous plot. Figures 3(b), (c), 
and (d) give, respectively, the plots corresponding to variations 4, 2, and 3 of the Mazzotta 
et al. (1998) DR rate coefficients. Variations 4, 5, 7, and 9 are very similar as are variations 
1, 2, 6 and 8. The curves due to Si XII 499.41 A have been omitted from these plots due to 
suspicions about the instrument calibration. In all panels the degree of overlap between the 
various intersections of the emission measure loci has improved. The improvement is best for 
variation 3, followed by variation 2. Variation 4 is the least successful. Thus the important 
conclusion that we are observing an isothermal plasma is, if anything, strengthened by our 
considerations of these plausible variations in the DR rate coefficients. 



5. Conclusions 

We have described on the basis of atomic physics theory what uncertainties may be 
present in the DR rate coefficients commonly used in ionization balance calculations. Of 
course this situation will only be definitely improved by further calculations and experiments, 
but given the number of individual rate coefficients required for astrophysical modeling 
purposes, this is an enormous amount of work. To help prioritize the needed atomic data, 
and to provide further motivation, we have investigated how the uncertainties in the DR 
rate coefficients translate into uncertainties in inferred FIP factors. We find that, depending 
on the specific line ratio, inferred FIP factors can be a factor of 5 smaller or 1.6 times larger 
that the FIP factor inferred using the unsealed DR rate coefficients of Mazzotta et al. (1998). 

Taking the unconstrained data in Table 2 and using the ratio of the maximum over 
minimum value of S for each line ratio, we can prioritize the need for DR measurements and 
calculations of the various isoelectronic sequences studied here. Listing the various line ratios 
in decreasing order of S max /S min yields: Si IX/S IX, Mg IX/S IX, Si X/S X, Mg IX/S X, 
Si IX/S X, Mg VI/Ne VI, and Mg VII/Ne VII. Based on the order and frequency with which 
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the different isoelectronic sequences occur in this list, we prioritize our list of needed rate 
coefficients for DR onto specific isoelectronic sequences, in decreasing order of importance, 
as follows: O-like, C-like, Be-like, N-like, B-like. To this we append F-like, Li-like, He-like, 
and Ne-like based on the range of the DR scale factors given in Table 4. 

We have arrived at our list through admittedly subjective means, but we believe that 
in this case the ends justify the means. Faced with the current degree of uncertainty in the 
dielectronic recombination rate coefficients relevant to astrophysical plasmas, the number of 
required calculations and measurements is a daunting task. Our aim in this work has been 
to prioritize this work, and to point out those ions where the uncertainties in the DR rate 
coefficients has the most impact on the analysis of astrophysical UV and X-ray spectra. 
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Temperature (K) 



Fig. 1. — Intensity curves for the emissivity (photons per second per unit electron density 
per ion) times ionization fraction of the Mg VI 1190.09 and Ne VI 558.59 A lines. Both 
curves have been normalized to their maximum values. The correspondence between the 
two is almost exact for T < 5 x 10 5 K. At higher temperatures, the Mg VI line is stronger. 
A steeply rising emission measure may weight emission from this temperature region enough 
to produce large apparent Mg/Ne abundance ratios, if this difference between the two curves 
is not accounted for during data analysis. 
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Fig. 2. — Emission measure loci for the ions Si VII- Si XII determined from SUMER ob- 
servations of an equatorial streamer using the ionization balance calculations of Mazzotta 
et al. (1998). The curves intersect at a temperature close to logT = 6.15, suggesting the 
conclusion that the plasma is isothermal. The spectral lines considered for each ion are: 
Si VII 1049.22 A; Si VIII 944.38, 949.22, and 1445.75 A; Si IX 950.14 and 694.70 A; Si X 
638.94 A; Si XI 580.91 A; and Si XII 499.41 A. These give three curves for Si VIII, two for 
Si IX and one for each remaining ion. 
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Fig. 3. — Emission measure (EM) loci from Figure 2 for the temperature region 6.0 < 
logT < 6.3. Figure 3(a) shows the results using the Mazzotta et al. (1998) ionization balance. 
Figures 3(b), 3(c), and 3(d) show, respectively, the plots corresponding to variations 4, 2, 
and 3 of the Mazzotta et al. ionization balance. Variations 4, 5, 7, and 9 are very similar as 
are variations 1, 2, 6 and 8. As can be seen from the Figure, variation 3 appears to give the 
best improvement to the inferred (and expected) isothermal nature of the observed coronal 
plasma, followed by variation 2, with variation 4 being the least successful. The curves due 
to Si XII 499.41 A have been omitted from these plots due to suspicions about the instrument 
calibration. 
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Table 1. Scale factors for the recommended DR rate coefficients of Mazzotta et al. (1998) 

onto selected ions of Ne, Mg, Si, and S. 



Isoelectronic 
Sequence 


Ions 


DR Scale Factors 


He 


Ne IX, Si XIII 


- 0.8 


1.0 1.2 


Li 


Ne VIII, Mg X, Si XII 


- 0.8 


1.0 1.2 


Be 


Ne VII, Mg IX, Si XI 




1.0 1.6 


B 


Ne VI, Mg VIII, Si X 




1.0 1.7 


C 


Ne V, Mg VII, Si IX, S XI 


.31 .62 


1.0 - 


N 


Mg VI, Si VIII, S X 


- .45 


1.0 - 





Mg V, Si VII, S IX 




1.0 3.6 


F 


Si VI, S VIII 




1.0 4.7 


Ne 


Si V 


- .85 


1.0 - 



Table 2. Range of scale factor S for all 1728 sets of DR variations (unconstrained) and for 

the reduced set of 274 variations (constrained). 



Line ratio 




Scale 


Factor S 




Unconstrained 


Constrained 


Minimum 


Maximum 


Minimum 


Maximum 


Mg VI/Ne VI 


0.60 


1.11 


0.62 


1.11 


Mg VII/Ne VII 


0.67 


1.22 


0.75 


1.22 


Mg IX/S IX 


0.33 


1.29 


0.41 


1.29 


Mg IX/S X 


0.51 


1.64 


1.00 


1.64 


Si ix/s ix 


0.20 


1.01 


0.24 


1.00 


Si ix/s x 


0.36 


1.14 


0.43 


1.14 


Si x/s x 


0.43 


1.60 


0.71 


1.60 
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Table 3. Observed Coronal FIP Fractionations (from Laming et al. 1999) 



Line ratio 


position l a 


position 2 


position 3 


position 4 


Mg VI 1190.09/Ne VI 558.59 




7.4 ± 


1.4 


4.5 ±0.8 


3.5 ±0.8 


Mg VII 868.11/Ne VII 895.17 


2.0 ±0.4 


2.6 ± 


0.5 


4.1 ±0.7 


4.1 ±0.7 


Mg IX 749.55/S IX 871.71 


3.6 ±0.8 


3.5 ± 


0.4 


3.2 ±0.1 


3.1 ±0.1 


Mg IX 749.55/S X 776.37 


3.0 ± 1.1 


3.2 ± 


1.0 


2.6 ±0.8 


3.0 ±0.9 


Si IX 950.14/S IX 871.71 


4.8 ±2.4 


3.6 ± 


0.5 


2.7 ±0.1 


2.8 ±0.1 


Si IX 950.14/S X 776.37 


4.6 ± 1.7 


3.0 ± 


0.9 


2.0 ±0.7 


2.4 ±0.8 



a Position 1 corresponds to a slit position covering 14" in radial distance over the 
solar limb. Positions 2, 3, and 4 are successively 14" further out in radial distance 
from the solar limb. 



Table 4. Scale factors for initial ions by which the recommended DR rate coefficients of 
Mazzotta et al. (1998) have been multiplied by for the 9 variations selected in a "trial and 

error" fashion. 



Variation Isoelectronic Sequence 



Number 


He 


Li 


Be 


B 


C 


N 





F 


Ne 


1 


1.0 


0.8 


1.0 


1.7 


0.62 


1.0 


1.0 


1.0 


1.0 


2 


1.0 


0.8 


1.0 


1.7 


0.62 


1.0 


1.0 


4.7 


1.0 


3 


1.0 


0.8 


1.6 


1.7 


0.62 


1.0 


1.0 


4.7 


1.0 


4 


1.0 


1.0 


1.0 


1.0 


0.62 


1.0 


1.0 


1.0 


1.0 


5 


1.0 


1.0 


1.0 


1.0 


0.62 


1.0 


1.0 


4.7 


1.0 


6 


1.2 


0.8 


1.0 


1.7 


0.62 


1.0 


1.0 


1.0 


1.0 


7 


1.2 


1.0 


1.0 


1.0 


0.62 


1.0 


1.0 


1.0 


1.0 


8 


1.2 


1.0 


1.0 


1.7 


0.62 


1.0 


1.0 


1.0 


1.0 


9 


1.2 


1.2 


1.0 


1.0 


0.62 


1.0 


1.0 


1.0 


1.0 



